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Abstract. We present a comparison of the results of the ASCA survey of LMXB with the Inogamov & Sunyaev 

theory of accretion flow spreading on the surface of neutron stars. The ASCA survey of LMXB of Church & 

1^-^ . Balucihska-Church (2001) revealed a systematic variation of the luminosity of blackbody emission from the neutron 

' star spanning 3 decades in total X-ray luminosity suggesting that the level of blackbody emission is controlled 

^SJ , by the physics of the inner disk/stellar interface, which we can hope to understand. Two types of explanation 

■ exist: firstly that there is radial flow between the inner disk and star at all vertical positions above the orbital 

plane so that the height of the disk directly determines the area of star emitting. Secondly, the height of the 

emitting region on the star is not directly related to the disk properties but depends on the mass accretion rate 

as suggested by Inogamov & Sunyaev (1999) in their theory of accretion flow spreading on the stellar surface. We 

O [ find that the survey results for the emitting area agree with this theory at the lowest luminosities. However, for 

higher luminosities, the blackbody emission is stronger than predicted by spreading theory suggesting that the 

^ ' emitting area is controlled by radial flow between disk and star. 

■ 
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d 1. Introduction the blackbody component in LMXB spectra. Opinion has 

„ , ,r • • ■ T n T been divided between location on the surface or bound- 

ihe nature oi the A-ray emission m Low Mass X-ray , j_ ^^ ■ _li • j_- t ^ ^ 

. /TUT^^nN 1 r 1 • i fl-ry layer at the star, or m the inner accretion disk, and 

Binaries LMXB has tor many years been controversial. , i • , i- ^ ^ mi i 

^ , , ' n -1 this controversy continues to recent times, ihere are also 

It IS clear that the accretion flow deposits its energy partly i-ij-o- i. ji-i-v 

. . 1 • physical dinerences between models involving X-ray emis- 

m the accretion disk resulting m heating and X-ray emis- • • i i i i • i • ^ n 

, , r ri sion m the boundary layer and those involving the stellar 

sion, and partly m a boundary layer at the suriace ot the r ^ i i \ t r • i ^ 

' , ,1 r • in 1 1 r , surlace (see below). In the former case, it is assumed that 

star or on the stellar surface itseli, but the nature ot the , , • • j- i j- ^ j. j. 

. . 1 1 1 1 1 the angular momentum m the inner disk adjusts to that 

emission has not been agreed. Lirstly, there has been the r,, , . , , , -^uj-i u-uu 

' or the star m a boundary layer m the disk which becomes 

question or how many continuum components there are. , , r v oni • j, , i i , • , j 

™ r 1 1 1 . , XI, ^^rr^ i i i • i the source ot X-rays. ihirdly, the location, geometry and 

1 he spectra or the less bright LMXB are clearly dominated -r,! - i i , -i 

. . , . , , 1 r- 1 V, 1-1 size ot the Comptomzmg region have been controversial, 

by Comptonization which led to the use OI the Generalized -.^ ... , r iz-i ^-i- ii- 

T 1 , n-r • 1 n • Many theoreticians have favoured Comptonization taking 
i hernial model L -e ^/^co httmg ot bjxosat spec- 



tra (White et al. 1988). Only in the case of brighter sources 

did spectral fitting demand that a blackbody component , , , , , , , , , , , . , , , . , , 

, , , ^ ^, . , , ^ , ,,r,cr,\ 1 1 1 to these problems although this has not been widely ap- 

was added to this. Mitsuda et al. (1989) proposed a model • ; , mi ■ r ^ 



place in an inner, hot region close to the neutron star. 
However, the dipping class of LMXB provides answers 



predated. This group of ~10 sources exhibit decreases 
in X-ray intensity at the orbital period due to absorp- 
tion in the bulge in the outer disk where the accretion 
flow from the companion star impacts. Emission models 
are more strongly constrained in these sources, as mod- 
els have to fit not only the spectrum of the quiescent 
Send offprint requests to: M. J. Church source as in other types of LMXB, but also several lev- 



consisting of multi-colour disk (MCD) blackbody emis- 
sion from the inner accretion disk plus blackbody emis- 
sion from the neutron star, Comptonized in a region local 
to the star, and this two-component model is able to fit 
many sources well. Secondly, there is the identification of 
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els of dipping. Spectral evolution is not simple, and can- 
not be fitted by absorbed one-component models. It has 
been shown by analysis of data from many observations 
with Exosat, ASCA, BeppoSAX and Rossi-XTE that these 
sources are well-described by a model consisting of point- 
like blackbody emission plus Comptonized emission from 
an extended accretion disk corona (ADC) (e.g. Church et 
al. 1997, 1998; B alucinska- Church et al. 1999, 2000, 2001; 
Smale et al. 2001). The blackbody is absorbed rapidly 
in dipping showing that it is point-like, and has tempera- 
tures between 0.9 and ~2 keV. The Comptonized emission 
is shown by spectral analysis to be removed gradually, 
and the complex spectral evolution is well-described by 
progressive covering of this component in dipping by the 
absorber, the covering fraction rising smoothly in many 
cases from zero to unity. This, of course, requires the 
Comptonizing region to be extended allowing the size to 
be estimated. Moreover, measurement of dip ingress times 
directly provides the size of the region, and values of the 
radius of this region for several sources typically of 50,000 
km have been determined (Church 2001). It is clear that 
a very extended Comptonizing region with high electron 
temperature can only be identified with an ADC. The dip- 
ping sources also indicate that this must be a thin corona, 
since if an ADC of this large radius was spherical, it is un- 
likely that the absorber on the outer disk would be able 
to overlap the source region completely to produce the 
100% deep dipping observed in many sources (Smale et 
al. 2001). It is likely that the ADC in all dipping LMXB 
is very extended, which can probably be related to the 
strong effects of the neutron star forming the ADC by 
evaporation of the accretion disk. Moreover, it is not ex- 
pected that the dipping sources are atypical, and so an 
extended ADC is expected in all LMXB. The large mea- 
sured size of the ADC has several significant consequences 
as discussed below. 

Firstly, these measurements of the ADC radius tadc 
allow us to rule out Comptonization models that involve a 
localized region in the neighbourhood of the neutron star. 
Secondly, the ADC size measurements have an impor- 
tant consequence for blackbody emission observed from 
LMXB. It is expected that a substantial fraction of the 
total available energy will be dissipated initially in the ac- 
cretion disk, implying a large disk blackbody component. 
Howcivcr, because the ADC is vciry extcmded, all of the hot 
inner regions of the disk will be covered by ADC, and con- 
sequently it is expected, given the high optical depth of 
the ADC (Church 2001), that all disk blackbody emission 
will become Comptonized in the ADC so that no black- 
body emission from the disk itself should be visible. It is 
also clear on the basis of the thin, flat disk and thin, flat 
ADC geometry, that thermal emission from the disk must 
be the major source of seed photons for Comptonization 
in the ADC. Thus, blackbody emission clearly identified 
in the spectra of many LMXB (e.g. Church & Balucinska- 
Church 2001, see below) must be identified with the neu- 
tron star/boundary layer without any contribution from 
the disk being expected. This contrasts with a recent ten- 



dency to fit a disk blackbody plus comptt Comptonization 
model routinely to BeppoSAX data (e.g. Guainazzi et al. 
1998; SidoU et al. 2001; Oosterbroek et al. 2001). This is 
thus inconsistent with the implications of the large mea- 
sured sizes of Accretion Disk Coronae. 

The large ADC size is also relevant to the comptt 
Comptonization model used in the above fitting. It has 
been claimed that this is preferable at low energies 
(< 1 keV) to a cut-off power law representation of 
Comptonization as the latter implies that the spectral flux 
density continues to increase with decreasing photon en- 
ergy, even for very low energies below 0.1 keV. Clearly, at 
very low energies the flux must decrease because of the 
lack of seed photons. The comptt model due to Titarchuk 
(1994) calculates a Comptonized spectrum assuming a sea 
of soft thermal photons and so avoids an increase of flux 
at the lowest energies. It uses the Wien approximation to 
the seed photon spectrum, valid for hv » kT. However, 
usage of this model in which values of kT from spectral fit- 
ting as high as 1 keV and 2 keV are derived (e.g. Guainazzi 
et al. 1998) will contravene the assumptions of the model, 
and will overestimate the seed photon spectrum, and thus 
the Comptonized spectrum by as much as 100% at 1 keV. 
Thus a model with kT = 2 keV will only be valid above 
10 keV. Regarding the availibility of seed photons, we can 
calculate whether there is any shortage of seed photons be- 
low 1 keV as follows. The large size of the ADC referred to 
above ensures that seed photons produced in the accretion 
disk out to a radius of ^^50,000 km will be Comptonized 
as the disk is covered by ADC out to this radius. We have 
calculated the multi-colour disk blackbody spectrum of 
this part of the disk assuming a thin disk temperature 
profile T{r) by integration in radius (for appropriate val- 
ues of mass accretion rate) . It is clearly dominated by the 
soft photons from outer radii. For a luminosity of 10^® erg 
s^^, the photon spectral flux density forms a broad peak 
between 0.001-0.01 keV, and for a luminosity of 10'^* erg 
s~^, the flux peaks between 0.01 0.1 keV, showing that 
there is no deficiency of low energy photons. Thus, the 
cut-off power law is a perfectly good description of the 
Comptonized spectrum to energies at least as low as 0.1 
keV, whereas application of the comptt model without re- 
stricting kT for the seed photons to suitably low values 
will be invalid. 

Finally, the extended size of the ADC has implications 
for reflection in LMXB. Because all of the inner disk is 
covered by ADC, illumination of the disk by the neutron 
star will not be possible, and "hot" reflection of this source 
from the ADC will take place. This will not however, have 
the absorption features of a reflection component from the 
accretion disk. Moreover, although the ADC itself can il- 
luminate the disk and lead to reflection, the high values of 
optical depth obtained for the ADC (Church 2001) suggest 
that this component will be reprocessed in the ADC and 
so not observed. This may explain the lack of detections 
of reflected components in LMXB. In a black hole binary 
such as CygX-1, the ADC radius is much smaller (Church 
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2001), and so a reflection component can be expected, as 
observed. 

1.1. The ASCA survey 

The above two-component model used to explain dipping 
LMXB spectra has been applied recently to a survey of all 
classes of LMXB, including the Atoll and Z-track sources 
(Church fc Balucinska-Church 2001), and found to fit all 
sources well. As part of this survey, it was shown that 
the blackbody emission cannot originate in the accretion 
disk, since fitting a two-component model consisting of 
multi-colour disk blackbody plus a Coniptonization term 
gave in many cases a value for the inner disk radius n 
more than 10 times smaller than the radius of the neu- 
tron star assumed to be ^ 10 km. Thus the blackbody 
emission has to be from the neutron star and these re- 
sults support the idea that disk blackbody emission is fully 
Comptonizcd in the ADC. The main result of the survey 
was the systematic behaviour that it revealed in the neu- 
tron star blackbody emission as reproduced in Fig. 1. In 
brighter sources, the 1 30 kcV blackbody luminosity Lbb 
approached 50% of the total luminosity (the dotted line) 
as expected in the Newtonian approximation. In fainter 
sources, -Lbb fell to much smaller fractions of the total. It 
was demonstrated that, as the total luminosity increased, 
the blackbody emitting area increased by a factor of 5400 
across the sample, whereas the blackbody temperature 
kTsB changed (in fact, decreased) by 40% (i.e. a factor 
of 0.60), corresponding to a change by a factor of 8 in T^. 
Thus, the area is the most important factor in determining 
the level of blackbody emission. 
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Fig. 1. Variation of blackbody luminosity in the band 1- 
30 kcV with total luminosity in the same band from the 
ASCA survey of LMXB of Church & Balucinska-Church 
(2001) 



1.2. Possible modification of the blackbody 

The survey assumed that there was simple blackbody 
emission from the surface of the neutron star, and we 
next discuss this assumption. Firstly, modification of the 
blackbody was proposed by Rutledge ct al. (1999) for 
transiently-accreting LMXB in quiescence, having kT typ- 
ically 0.2 keV and luminosity L ~ 10^^ erg s~^. The mod- 
ification arises from higher energy photons escaping from 
greater depths in the atmosphere. This is clearly inappli- 
cable to LMXB not in quiescence emitting at luminosities 
of 10^^-10'^* erg s~^ (Rutledge priv comm). 

Secondly, although in the study of X-ray bursts there 
is observational evidence that there is sometimes mod- 
ification of the blackbody by electron scattering in the 
neutron star atmosphere during bursts, there is no such 
evidence in the case of non-burst emission. Fitting a sim- 
ple blackbody in a case where modification did take place 
would lead to errors in kT and the emitting area. There 
are two pieces of evidence for modification in bursts (see 
the detailed discussion in Balucinska-Church et al. 2001). 
Firstly, values of kT of ~3 keV have been obtained for a 
fraction of bursts from fitting a simple blackbody to the 
peak, which implies super-Eddington emission. Secondly, 
an increase in blackbody radius it!BB while the blackbody 
temperature decreased during burst decay was found in 
the source XB 1636-536 (Inoue et al. 1984; Sztajno et al. 

1985) . The increase of i?BB was thought to be unreal and a 
consequence of fitting a simple blackbody when the emis- 
sion was modified blackbody. However, whether this effect 
is observed depends on whether spectral fitting allows for 
blackbody emission by the non-bursting source (Sztajno 
et al. 1986), and so is not strong evidence for modifica- 
tion. More recently, Kuulkers et al. (2002) have presented 
evidence that during the cooling phase of burst decay in 
GX 17-1-2, the blackbody is not modified. 

In fact, whether modification takes place depends crit- 
ically on the electron density rie which determines whether 
absorption processes or electron scattering dominate the 
opacity. Various theoretical descriptions of the neutron 
star atmosphere in bursting have been given, by Fujimoto 
et al. (1981), Paczyhski (1983), London ct al. (1984, 

1986) , Ebisuzaki et al. (1984), Ebisuzaki (1987) and Madej 
(1991). However, there are very large disparities between 
the electron density values in these works, rio varying be- 
tween 10^**- 10'^^ cm~'^. While an atmosphere with = 
10^* cm~^ is clearly dominated by electron scattering, a 
higher density atmosphere is not. Thus the theoretical pic- 
ture is unclear and we have to rely on observational ev- 
idence that some bursts only have a modified spectrum. 
For these, the effect of modification on measured values 
of kT depends on the ratio of colour temperature to ef- 
fective temperature Tcoi/Tes for which recent theoretical 
values are low, i.e. 1.1 in the work of Madej & Rozanska 
(2000) (Czerny, priv comm), compared with earlier values 
which were higher. Values can be derived using the London 
et al. model of ^^1.6 for fcTpff ~ 2 keV appropriate to a 
burst peak, ~1.3 for 1.5 keV, and ~1.1 for 1 keV (assum- 
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ing hydrogen composition). Thus even these factors imply 
relatively small corrections for non-burst emission. The ef- 
fect on emitting area depends on this ratio squared, so in 
terms of recent theory can reach only 20%. In non-burst 
emission, there is actually no evidence for modification at 
all and so it would be inappropriate to apply a correc- 
tion factor. On the contrary, the ASCA survey results in 
which the blackbody radius approaches 10 km in bright 
sources, i.e. consistent with the neutron star radius pro- 
vide evidence that there is no large modification. Thus it 
is unlikely that the results of the present paper would be 
affected. 
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Fig. 2. Variation of the height h of the blackbody emission 
region on the surface of the neutron star with the height of 
the inner, radiatively-supported accretion disk calculated 
at the radial position where Pr = 10 Pg 



In Fig. 2, the ASCA survey results are replotted in the 

form of half-height h of the omitting region on the neutron 
star assumed to be an equatorial strip, more specifically, 
a sphere intersected by two parallel planes, having emit- 
ting area An Rh whore R is the radius of the neutron 
star, h varied between 20 m (at the lowest luminosities) 
and the full radius of the star ~ 10 km. These heights 
are (for almost all of the data) substantially larger than 
the half- height of the inner accretion disk calculated using 
thin disk theory (Shakura & Sunyaev 1976) at some posi- 
tion close to the surface of the neutron star (e.g. l.OOli?). 
This thin disk half-height is ~20 m for a luminosity of 
10^^ erg and ~40 m for 10^^ erg s~^. However, the 
inner disk is expcx'ted to be radiativcly supported over a 
wide range of luminosities > 3 x 10^^ erg s~^ and conse- 
quently has much larger vertical extent. The radiatively- 
supported disk merges at some radius with the thin disk 
and the half- height H{rio) was calculated at the radius at 
which the radiation pressure was 10 times the gas pres- 
smc (Czerny & Elvis 1987). Remarkably, there was good 
agreement between h and H{r lo) for the majority of the 



survey sources, the agreement spanning 3 decades in each 
quantity as shown in Fig. 2 (Fig. 3, Church & Balucinska- 
Church 2001). This agreement suggests that the proper- 
ties of the inner disk determine in some way the blackbody 
emitting area. In the case of the bright source GX 5-1, the 
agreement is poor with h = 110 km and H(rio) = 43 km. 
The blackbody radius Rbb is 33 km showing that the data 
cannot be explained by emission from the neutron star 
alone, suggc;sting that the emitting region has expanded 
to a spherical cloud centred on the star. However, there 
is good agreement between the height of the cloud at 33 
km (= Rbb) plotted as a triangle in Fig. 2 and Hirio). 
A similar point is plotted for the blackbody radius of the 
source GX 13-1-1. There is thus some evidence that, for 
these points also, H determines h. In the present paper, 
we compare the survey results with possible mechanisms 
leading to the agreement of h = H, specifically i) radial 
flow between the accretion disk and neutron star, and ii) 
accretion flow spreading on the surface of the neutron star 
(Inogamov & Sunyaev 1999). 

The agreement between h and H suggests that the 
properties of the inner disk determine the blackbody emit- 
ting area in some way. However, there are at least two 
possible explanations of the agreement. Firstly, there may 
be radial flow across the gap between the inner disk and 
the star at all vertical positions above the orbital plane so 
that the disk height approximately sets the height of the 
emitting area on the star. This requires there to be a sta- 
ble, radiatively-supported inner accretion disk, although 
it has been suggested that the radiative disk is subject 
to instabilities. These may not however, disrupt the disk 
in a major way as their characteristic timescale is short 
(see Sect. 4). Secondly, there may be spreading of the ac- 
cretion flow on the neutron star surface as proposed by 
Inogamov & Sunyaev (1999). In this case, accretion flow 
between the disk and the star will be essentially in the 
orbital plane. In this mechanism discussed below, the ac- 
cretion flow spreads vertically over the neutron star, the 
height reached depending on the mass accretion rate and 
thus on the total luminosity. The emitting area depends 
on L, and the quantity H is simply a measure of L; thus 
the agreement of h with H represents the variaticm of 
emitting area with luminosity (and does not result from a 
direct mapping of the disk height onto the star). We next 
compare the survey results with this theory. 

2. Accretion flow on the surface of the neutron 
star 

In the often accepted standard theory of accretion flow 

interaction with neutron stars, it is assumed that the ad- 
justment in angular velocity between the accreting mate- 
rial with high Keplerian velocity and the slower rotation of 
the star takes place in a boundary layer located in the in- 
ner disk (Shakura & Sunyaev 1973; Lynden-Bell & Pringle 
1974; Pringle & Savonije 1979; Papaloizou & Stanley 1986; 
Popham et al. 1993; Bisnovatyi-Kogan 1994; Popham & 
Narayan 1995). This takes place via turbulent friction be- 
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tween differentially rotating layers in the disk. Thus X-ray 
emission would take place in this inner disk. 




0.1 1 10 



Arc length Re (km) 

Fig. 3. Variation of X-ray flux q emitted from the neutron 
star expressed as a fraction of the Eddington flux go with 

arc length R6 for 4 cases in which the source luminosity 
is (a) 1%, (b) 4%, (c) 20% and (d) 80% of the Eddington 
limit 



Inogamov & Sunyaev (1999) recently proposed a radi- 
cally different approach in which the adjustment in angu- 
lar velocity takes place on the surface of the star, mate- 
rial meeting the star in the equatorial plane and flowing 
vertically on the stellar surface away from the plane and 
forming a spreading layer. Viscous interaction with the 
underlying surface takes place while the angular velocity 
of the flow material differs from that of the star, load- 
ing to heating and X-ray emission which rises to a peak 
at a particular arc distance R9 measured on the surface 
of the star from the equator. The emitted flux then falls 
rapidly with R9 as energy is radiated. Thus, for a given 
luminosity, there is a band on the surface of the star where 
the X-ray emission takes place. A system of equations was 
set up for the spreading layer considered essentially as a 
stellar atmosphere, consisting of hydrostatic equilibrium, 
radiative equilibrium, an equation of state containing gas 
and radiation pressure and the Stefan-Boltzmann law. The 
opacity was assumed to be dominated by Thomson scat- 
tering. Transformation and solution of these equations was 
carried out for a number of values of accretion rate, i.e. 
neutron star spreading layer luminosity. In particular, re- 
sults were presented for the variation of X-ray flux q de- 
fined as the energy emitted per unit area per second, with 
arc distance. Fig. 3 reproduces the results of Inogamov 
& Sunyaev for the distribution of X-ray flux with R 9 for 
4 values of the neutron star luminosity equal to 1%, 4%, 



Table 1. Emitting region on the surface of the neutron 
star calculated using the theory of Inogamov & Sunyaev 
(1999) expressed in terms of the range oi R9 and also in 
terms of half-heights h derived for these for 4 luminosities 
of the spreading layer Lsl 





ARO 


h 


ARO 


h 






FWHM 




7-0.99 
^0.5 


% 


km 


km 


km 


km 


1 


0.12 


0.12 


0.09 


0.09 


4 


0.44 


0.44 


0.32 


0.32 


20 


2.47 


2.45 


1.51 


1.50 


80 


15.81 


11.62 


9.14 


8.28 



20% and 80% of the Eddington limit LEdd- The flux is 
given relative to the Eddington flux qo, where 

LEdd = inR^o = 1-26 x lO^^ -^erg s'^, 

Mq 

and 

_ GMnipC 
Qo — — 52 

it'' (TT 

where M is the mass of the neutron star, rup is the proton 
mass and ctt is the Thomson cross section. i^Edd is, of 
course, independent of the radius assumed for the neutron 
star. Fig. 3 shows q/qo integrated over photon energy as a 
function of arc length assuming qo appropriate to a 1.4Mq 
object so that L^dd is 1-76 x 10^^ erg s~^. For the 4 cases, 
the neutron star luminosities are 1.76 x 10^^, 7.04 x 10^^, 
3.52 X 1037 and 1.67 x lO^^ erg s'^ 

Clearly, at low luminosities the X-ray emission peaks 
in a narrow strip relatively close to the orbital plane, but 
extending down to = 0. At high luminosities the hot re- 
gion expands to fill the star. Inogamov & Sunyaev showed 
that the calculated fluxes were not sensitive to the mag- 
nitude of the viscosity parameter, as might be expected. 
X-ray observations are not capable, of course, of revealing 
the values R9i and R92 between which most of the X-ray 
emission is concentrated assuming the model to be cor- 
rect, but they do provide values of the emitting area from 
fitting a simple blackbody. From the measured luminos- 
ity and temperature, an area can be derived and compared 
with the emitting area of the neutron star in the Inogamov 
& Sunyaev theory, as described in the next section. 

3. Results 

Fig. 3 shows the variation of flux with polar angle q{R9) 
expressed as a fraction of the Eddington flux qo in the 
theory of Inogamov & Sunyaev which we wish to compare 
with the results of the ASCA survey of LMXB. Although 
Inogamov & Sunyaev use an electron scattering dominated 
opacity for the surface layers of the neutron star, they in 
fact, implement their model assuming simple blackbody 
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emission. As discussed in Sect. 1, whether modification of 
the blackbody takes place depends on the electron density, 
and no evidence for modification in non-burst emission 
is known to exist. Thus we simply obtain the theoretical 
sizes of the hot regions on the neutron star from Fig. 3 
and compare these with blackbody areas from the survey. 

Wc adopt two ways of defining the emitting area from 
the theoretical variation q{R9); comparison of these will 
indicate the uncertainty in the area. Firstly, we use the 
region on the neutron star corresponding to the FWHM 
of the flux q(R0), i.e. the range between points where q 
= 9max/2. Secondly, we used a definition based on the in- 
tensity I oc f q{6 ) cosd dd integrated over the surface of 
the sphere allowing for the angular dependence of the area 
2-jtR'^ cos6 d6. The intensity I{6) increases asymptotically 
to a constant value corresponding to the maximum value 
oi R9 at which there is X-ray emission. The range oi R9 
was then obtained between I = 0.50 /max and 0.99 /max- 
For each method, we calculated the half-height h of an 
equatorial strip having the same area, from the range of 
R6 values using h = RsinO with 12 km for the radius 
of the neutron star as assumed by Inogamov & Sunyaev. 
These data are shown in Table 1. The case of the highest 
luminosity point is interesting. The FWHM in this case 
gives a value of A{R9) of 15.8 km which corresponds to 
A9 = 75.5° if we assume R ^ 12 km, and 90.5° for R = 
10 km. Thus, in the latter case, the whole star is clearly 
emitting and the value of h equal to Rsm9 becomes 10 km. 
Alternatively, if we assume that 15.8 km is one quarter of 
the circumference of the star, a radius of 10.1 km for the 
neutron star is obtained. 
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Fig. 4. Comparison of 4 theoretical points for the ver- 
tical extent of the emitting region on the neutron star 
with ff(rio) (solid curve), the height of the disk where Pr 
= lOpg. Open squares and open circles are uncorrected 
(see text); the squares show the height equivalent to 99% 
of intensity, the circles show the height equivalent to the 
FWHM of the flux distribution. Filled squares and fllled 
circles are at corrected values of the total luminosity 



In Fig. 4 we compare values of h calculated in these 
two ways with H{rio), the height of the inner disk, as a 
function of total luminosity. This function was shown to 
agree well with the ASCA survey results by Church & 
Balucihska-Church (2001), and was calculated assuming 
a neutron star radius of 10 km. Open squares show the 
99% intensity points, and open circles show the FWHM 
points. It can be seen that there is less than a factor of two 
difference! between the theoretical points of Inogamov & 
Sunyaev for the two methods of area calculation. However, 
the contribution of the neutron star blackbody is only 
part of the total luminosity in X-ray binaries and so the 
theoretical points should be corrected. The other major 
emission component is that arising originally as thermal 
emission in the accretion disk, which is then Comptonized 
in the ADC (e.g. Church 2001; Balucihska-Church et al. 
2001) making a dominant contribution in most LMXB. 
As Inogamov & Sunyaev (1999) assumed the Newtonian 
approximation in which the energy available at the inner 
edge of the accretion disk is one half of the total energy, 
one possible method of correction would be to plot the 
points at Lxot — 2-Lbb, where Lbb is the luminosity of 
the spreading layer. The points corrected in this way are 
shown as filled circles and squares. The differing values of 
neutron star radius used by Inogamov & Sunyaev (12 km) 
and in the calculation of H{rio) (10 km) will have only 
10% effect on the comparison. 




1 10 100 

Lr„t (10^ erg s"') 



Fig. 5. Same as Fig. 4 with data points from the LMXB 
survey added for luminosities in the band 1-30 keV. The 
set of points labelled A are corrected as in Fig. 4 The 
points labelled B arc corrected in luminosity based on 
observational information from the LMXB survey. Also 
shown is the variation of the thin disk half-height cal- 
culated close to the stellar surface (at 1.001 i?) with L. 
The theoretical height of the radiatively-supported disk 
(curved line) approaches the thin disk line at low L where 
radiative support fails 
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In Fig. 5 we add the survey data points to the com- 
parison. The theoretical points are shown as before, un- 
corrected and also corrected by a factor of two shift in 
i'Tot • We now consider the extent of the required correc- 
tion further. Inogamov & Sunyaev (1999) assume that the 
spin of the neutron star is essentially zero, and SibgatuUin 
& Sunyaev (1998) argue that in this case, Lbb = 2.8-iTot- 
This implies that Lxot = 1.34 -Lbb, so that we should cor- 
rect the theoretical points in Figs. 3 and 4 by moving them 
to higher values of Lxot by 34%. However, this argument 
predicts that the spectra of LMXB in general would be 
strongly dominated by the thermal emission of the star, 
whereas, in fact, in many sources the Comptonized emis- 
sion strongly dominates, and work on the dipping sources 
has proven that this originates in an extended ADC. It 
has been suggested, for example by Inogamov & Sunyaev 
(1999) that modification of the neutron star emission in 
the stellar atmosphere will take place; however this will 
not lead to a power law spectrum as observed, but to a 
spectrum very similar to a simple blackbody (Madej 1991) 
and the observed dominance of Comptonized emission in 
the spectrum could only be avoided by assuming that most 
LMXB have a neutron star spinning sufficiently fast that 
the emission due to viscous dissipation in the surface is 
markedly reduced. As an alternative to seeking a theoreti- 
cal relation between Lxot and Lrb, we can use an approx- 
imate observational relationship provided by the ASCA 
survey. Fig. 1 shows that there is not a simple relation- 
ship between Lbb and Lxot , but rather that the blackbody 
luminosities occupy a broad band of the two-dimensional 
space. To make an approximate correction, we use the cen- 
tre of this band. This requires a substantial movement of 
the theoretical points in Fig. 5 to higher values of ixot as 
shown by the filled squares labelled B. 

However, we can avoid the difficulty of the correction 
by comparing the theoretical points with the survey data 
replotted in terms of Lbb instead of Lxot as shown in Fig. 
6. A curve showing the theoretical height of the accretion 
disk is not shown as this is a function of Lxot- It can be 
seen that the survey points tend to converge with the the- 
ory at the lowest luminosity, but are substantially larger 
than the theoretical points at higher luminosities. 

4. Discussion 

The main result of the ASCA LMXB spectral survey was 
the agreement oi h = H, for sources in the quiescent state 
(i.e non- flaring) . However, it is not clear which of two ex- 
planations of this is correct: that of radial flow between 
inner disk edge and star, in which the height of the inner 
radiativcly supported disk directly determines h, or the 
theory of Inogamov & Sunyaev (1999) in which spreading 
of material on the neutron star surface to height h depends 
on the mass accretion rate of which iJ is a measure. In 
the case of the first mechanism, the inner disk may have 
the steep radial profile implied by standard theory close 
to the neutron star, but radial flow takes place across the 
gap between inner disk and star thus defining the emitting 
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Fig. 6. The 4 theoretical points plotted as a fuction of 
blackbody luminosity compared with data from the survey 
in the band 1-30 keV 



region on the star of similar height to the highest point of 
the inner disk. Although there; has been extensive theoret- 
ical investigation of advective flow in black hole systems 
(e.g. Abramowicz et al. 1996), very little work has been 
carried out for neutron star binaries, and it is not clear 
whether the sonic point lies within the stellar surface or 
not. However, Popham fc Sunyaev (2001) show that the 
radial velocity increases by two orders of magnitude within 
the disk. Two- or three-dimensional modelling would be 
required to decide whether radial flow across the gap can 
take place. A possible objection to the first mechanism is 
that the radiatively supported inner disk may not be sta- 
ble, i.e. implying that it cannot exist. The stability of the 
inner disk in a Galactic black hole system GRS 1915+109 
has been considered theoretically by Janiuk et al. (2000) 
who showed that oscillations are expected leading to peri- 
odic variations in the X-ray luminosity of the system sim- 
ilar to those observed. However, it was not found that this 
leads to major disruption of the disk. Since the timescale 
for variability depends on mass (Czcrny, priv. comm.), we 
have applied this model to LMXB by scaling the stellar 
masses, and find a timescale of ~ 30 s. Variability on this 
timescale is not obviously observed, and would not in any 
case prove that major disruption of the disk takes place. 

In the present work, we have compared the survey re- 
sults with the theory of Inogamov & Sunyaev (1999) in 
two ways. Firstly, we made an approximate comparison as 
shown in Figs. 4 and 5. When a more appropriate compari- 
son is made by comparing with Lbb, not Lxot, as shown in 
Fig. 6, the survey results and theory converge at the low- 
est luminosities, but at higher luminosities, the data, i.e. 
the emitting areas, exceed the theory by a factor that we 
estimate as between four and eight. However, the methods 
used of obtaining the emitting area from theory underes- 
timates the area, for example by using the FWHM, by up 
to a factor of two. It consequently appears that the data 



8 



Church ot ah: Blackbody emission in LMXB 



exceed the theory by a factor of two to four. Thus at low 
luminosities, the emitting area on the neutron star may 
be determined by accretion spreading, although even here 
radial flow would be an alternative explanation. At higher 
luminosities, the area exceeds the spreading theory imply- 
ing that the area is determined by radial flow which leads 
to the observed agreement between h and H and which en- 
hances the contribution due to spreading. Although there 
is no observational evidence for any modification of the 
neutron star blackbody spectrum in the non-burst emis- 
sion of a LMXB, this would lead to an underestimation of 
the emitting area, and the disparity between the survey re- 
sults and spreading theory would be somewhat increased. 

Extension of the present accretion flow spreading 
theory which is essentially one-dimensional, to a two- 
dimensional form would be desirable. Further analysis of 
LMXB spectra is clearly needed to increase the size of 
the sample and to make more use of BeppoSAX data to 
constrain spectral fitting parameters better. However, the 
present results indicate that while the neutron star accre- 
tion flow spreading theory of Inogamov & Sunyaev may 
set a base level to the neutron star blackbody emission, 
as observed in the faintest sources, it underestimates the 
observed blackbody luminosity in brighter sources. This 
implies that radial flow dominates leading to the agree- 
ment between h and H. 
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